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\ In this talk I summarise recent results on the infrared exponents of SU^Ay-Yang-Mills theory. 

I discuss a self-consistent power law solution for the Dyson-Schwinger equations for general 
1 Pi-Greens functions in the infrared. The corresponding running coupling has a fixed point at 
zero momentum, which turns out to be universal and gauge invariant within a class of transverse 
Ph, gauges. When calculated on a torus the infrared exponents of the ghost and gluon propagators 

differ from the corresponding continuum solutions. They agree, however, well with results from 
lattice calculations. 
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The infrared behaviour of the Green's functions of SU(/V c )-Yang-Miils theory is related to 
confinement in several interesting ways. One example is the Kugo-Ojima confinement criterion, 
which is satisfied in Landau gauge if the dressing function of the ghost propagator is singular and 
the gluon propagator is finite or vanishes in the infrared [Jl], §]. A second example is provided 
by the colour Coulomb potential, which can be determined from the instantaneous part of the 
time-time component of the Coulomb gauge gluon propagator [||]. In this talk I summarise recent 
analytical results on the infrared exponents of the Green's functions of SU(A^ C )-Yang-Mills theory 
in Landau gauge [Q] and a class of gauges that interpolate between Landau and Coulomb gauge 
[Q]. Furthermore, I discuss the relation of these solutions with corresponding results from Dyson- 
Schwinger equations (DSEs) and lattice simulations on a compact manifold [Q]. 

1. Infrared exponents in the continuum: Landau gauge 

The basic idea to determine the infrared behaviour of one-particle-irreducible (1PI) Green's 
functions is to investigate their Dyson-Schwinger equations order by order in a skeleton expansion 
{i.e. a loop expansion using full propagators and vertices). It turns out that these Green's functions 
can only be infrared singular, if all external scales go to zero. Thus to determine the degree of 
possible singularities it is sufficient to investigate the DSEs in the presence of only one external 
scale p 1 <C Ag CD . In collaboration with R. Alkofer and F. Llanes-Estrada we showed by induction 
that the skeleton expansion is stable in the sense that higher orders generate terms that are equally 
or less singular than those of lower orders. Thus one can determine the infrared behaviour of these 
functions from their DSEs approximated to lowest order in the skeleton expansion. We found a 
self-consistent solution of the tower of DSEs for lPI-Green's functions, given by [Q] 

H^C/) ~ {p 2 ) (n - m)K . (1.1) 

Here F n,m (p 2 ) denotes the infrared leading dressing function of the lPI-Green's function with 
2n external ghost legs and m external gluon legs. The exponent k is known to be positive [0], 
independent of any truncation scheme of the DSEs. 

A special instance of the solution ( |1.1[ ) are the inverse ghost and gluon dressing functions 
r 10 (^) 2 ) = G 1 (p 2 ) and r°' 2 (p 2 ) = Z~ l {p 2 ), which are related to the ghost and gluon propagators 
via 

The corresponding power laws in the infrared are 

G(p 2 ) ~ (p 2 )-*, Z(p 2 ) ~ (p 2 ) 2K . (1.3) 

In this notation the Kugo-Ojima criterion translates to the condition K > 0.5. For a bare ghost- 
gluon vertex in the infrared, justified by lattice calculations []8|] and also in the DSE-approach [Q], 
one obtains K = (93 — \/l201)/98 0.595 [10, 11], which satisfies the criterion. This specific 



Z/V) = -2g>, D^)=(s, v -^2£!. (1-2) 



value of K is found to vary only slightly for a large class of possible dressings of the ghost-gluon- 
vertex [11]. Similar values have been determined from exact renormalisation group equations [|T"2|]. 



An interesting consequence of the solution ( |L~l| ) is the qualitative universality of the running 
coupling in the infrared. Renormalisation group invariant couplings can be defined from either of 
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the primitively divergent vertices of Yang-Mills-theory, i.e. from the ghost-gluon vertex (gh — gl), 
the three-gluon vertex (3g) or the four-gluon vertex (4g) via 

agh-gitf) = |_G 2 (/)Z(p 2 ) p ~° const/A^, (1.4) 
47T 

« 3 V) = jV'V)] 2 zV) const/A^ (1.5) 

47C 

« 4 V) = f^'W^V) const/A^. (1.6) 

Using the DSE-solution ( |1.1| ) it is easy to see that all three couplings approach a fixed point in the 



infrared. This fixed point can be explicitly calculated for the coupling (1.6). Employing a bare 
ghost-gluon vertex one obtains a gh ~ gl (0) « S.92/N c [JTTJ]. 

2. Running coupling between Landau and Coulomb gauge 

One might expect that gauge invariant features of Yang-Mills theory show up in a quantity like 
the running coupling. Thus it is interesting to determine the coupling together with the infrared 
exponents of the corresponding dressing functions also in other gauges. A particularly interesting 
class of (transverse) gauges is specified by the gauge condition adoAo + V • A = with the gauge 
parameter < a < 1. The values a = 0, 1 correspond to Coulomb- and Landau-gauge respectively, 
whereas non-integer values of a interpolate between these gauges. 

In collaboration with D. Zwanziger we showed that in these gauges two renormalisation group 
invariant running couplings can be defined from the ghost-gluon vertex [^p. In the Landau gauge 
limit, a = 1, these couplings unify to the expression (|L4|). In the Coulomb gauge limit, however, 
they constitute two invariant charges which are given by 

«coui(|k|) = |^Zj,f(|k|), (2.1) 

«/(|k|)^y^G 2 (|k|)^M, (2.2) 

where Z^ st is the instantaneous part of the time-time component of the gluon dressing function 
and Zg r t is the equal-time spatial gluon dressing function. The first coupling, G5 CO ui, is nothing else 
than the familiar colour Coulomb potential, which has been shown to be (almost) linearly rising 
[TJ, 14]. The second coupling, aj, has a fixed point in the infrared This can be seen from 



the Dyson-Schwinger equation of the ghost propagator in Coulomb gauge. In interpolating gauges 
with < a < 1, this fixed point can be calculated and we find the same value 05/ (0) ~ S.92/N c as 
in Landau gauge. Correspondingly, the ghost and gluon propagator obey an infrared power law 
with exponent K m 0.595 independent of the gauge parameter a. The observation of these gauge 
independent features may have the potential to bring us closer to a unified picture of confinement 
in Landau and Coulomb gauge. 

3. Numerical solutions for the propagators in the continuum and on a torus 

There are at least three caveats in comparing results from the continuum Dyson-Schwinger 
approach to those of lattice calculations. First, the quantitative aspects of the continuum solutions 
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p 2 [GeV 2 ] p 2 [GeV 2 ] 

Figure 1: The gluon propagator function Z(p 2 )/p 2 (left diagram) and the ghost dressing function G(p 2 ) 
(right diagram). Compared are the results of the continuum DSEs with two respective solutions of DSEs on 
tori of different volumes. 



depend on the details of the chosen truncation scheme, whereas the lattice calculations are ab initio. 
On the other hand, lattice calculations are carried out on a compact manifold, and therefore one 
has to deal with effects due to finite volume and lattice spacing. Finally, boundary conditions have 
to be chosen on the compact manifold. To quantify some of the latter two effects we formulated 
the Dyson-Schwinger equations on a torus without changing the truncation scheme [TJ|]. One 
would then expect to see differences to the continuum solution for small volumes, which disappear 
continuously when the volume is chosen larger and larger. However, for our solutions this is not 
the case. In fig. [I] we compare solutions for the gluon propagator Z(p 2 ) / p 2 (left diagram) and the 
ghost dressing function G(p 2 ) (right diagram) of the continuum with those from two different tori 
with volumes of V « 2500 fm 4 and V « 10 5 fm 4 . Note that the smaller volume corresponds to a 
typical one achieved in recent lattice simulations [[[(]], whereas the larger value exceeds the current 
possibilities of the lattice formulation by far. Both solutions on the torus agree very well with each 
other, whereas there is a clear difference to the continuum results. In terms of the infrared exponent 
K, defined in eq. (|L3|), one finds a value of K 0.5 for the gluon propagator on the torus, which 
should be compared to the corresponding continuum value K m 0.6. Similarly, whereas the ghost 
dressing function diverges in the continuum, we observe a finite ghost on the compact manifold. 
These differences persist, when the infrared behaviour on the torus is extrapolated to the continuum 

We compare these results to the ones of recent lattice simulations in fig. § and fig. [| The lattice 
simulations for the gluon propagator are all performed for SU(3) at different values of j8, different 
lattice sizes and different actions [16, 17, 18] (see also Jl9|]). They fall on top of each other. On 



the other hand there are slight differences between the SU(2) ghost dressing function of ref . [ 2C ] 



and the SU(3) result of ref. [|18|] (see also Q21|]). Interestingly, there is nice qualitative agreement 
between both, the lattice and the DSE-formulation on the compact manifold. Both approaches give 
a finite gluon propagator in the infrared. This is in marked difference to the numerical continuum 



result, which agrees with the analytic power law, eq.(1.3), with K 0.6. The situation is less clear 



for the ghost dressing function, though at least the SU(3) result of ref. [18] seems to follow more 
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Figure 2: The results for gluon and ghost from Dyson-Schwinger equations in the continuum and on the 
torus are compared with the lattice data of refs. flip] , |l7| , |l8[ 
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Figure 3: Results for the running coupling of the ghost-gluon vertex from DSEs in the continuum and on 
the torus compared to the lattice data of ref. [p"8||. 



or less the DSE-solution on the torus. For the running coupling 0C gh ~ gl (p 2 ) = ^G 2 (p 2 )Z(p 2 ) the 
numerical solution in the continuum reproduces the fixed point a gh ~ gI (0) as S.92/N c discussed 
below eq.(pTrJ). On the compact manifold, however, the coupling vanishes at small momenta 1 , a 
behaviour which seems to be at least counterintuitive. 

Our results from DSEs compared to lattice simulations of the ghost and gluon propagators 
show an interesting difference between the infrared behaviour found in the continuum and the 
one found on a compact manifold. Finite volume effects can hardly explain this result, since the 
difference persists for extremely large tori. The role of different boundary conditions in this respect 



is currently investigated (c.f. also ref. [24], these proceedings) 



A similar result has been found for the running coupling of the three-gluon vertex in refs. [§|,[j3|]. Note, however, 
that the conclusions drawn in [ p3| ] are contradicted by the results of refs. [Q, ^| p[p"i[]. 
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